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Abstract 

The leading excitations of the dilute Al model in regime 2 are considered 
using analytic arguments. The model can be identified with the integrable 1 2 
perturbation of the unitary minimal series Mll+i- It is demonstrated that the 
excitation spectrum of the transfer matrix satisfies the same functional equations 
in terms of elliptic functions as the exact S-matrices of the <j>i 2 perturbation 
do in terms of trigonometric functions. In particular, the bootstrap equation 
corresponding to a self-fusing process is recovered. For the special cases L = 3, 4, 6 
corresponding to the Ising model in a magnetic field, and the leading thermal 
perturbations of the tricritical Ising and three-state Potts model, as well as for the 
unrestricted model, L = oo, we relate the structure of the Bethe roots to the Lie 
algebras £"8,7,6 and D4 using Coxeter geometry. In these cases Coxeter geometry 
also allows for a single formula in generic Lie algebraic terms describing all four 
cases. For general L we calculate the spectral gaps associated with the leading 
excitation which allows us to compute universal amplitude ratios characteristic of 
the universality class. The ratios are of field theoretic importance as they enter 
the bulk vacuum expectation value of the energy momentum tensor associated 
with the corresponding integrable quantum field theories. 
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1 Introduction 



One of the most significant applications of two-dimensional integrable systems in field 
theory and statistical mechanics is to systems undergoing a continuous phase transition. 
Near the critical point the correlation length diverges rendering the system scale invari- 
ant. This leads to the well-known concept of universality classes, which describe the 
macroscopic scaling behavior of the thermodynamic quantities near the critical point. 
The universality classes can be identified with conformal field theories which encode the 
critical exponents in the spectrum of their primary fields. Besides the critical expo- 
nents, each class exhibits additional universal ratios of critical amplitudes, which might 
equally well serve as a characteristic, see e.g. and references therein. In two dimen- 
sions exact non-perturbative results are available on these critical amplitudes, provided 
the off-critical model can effectively be described by an integrable perturbation of the 
conformal field theory governing the critical point. There are various ways to proceed, 
such as conformal perturbation theory, the bootstrap approach involving the construc- 
tion of exact scattering matrices and the transfer matrix method. While we concentrate 
on the last of these in this article, it has intimate ties with the other techniques which 
we will point out in the course of our discussion. 

We consider the integrable off-critical dilute Al lattice model @ which belongs to 
the class of interaction-round-a-face (IRF) models with adjacency matrix related to the 
simple Lie algebra Al. The epithet "dilute" points out that neighbouring lattice sites are 
allowed to have the same height value in addition to those permitted in the conventional 
IRF models, where neighbouring heights take adjacent values in the Dynkin diagram. 
The dilute Al model possesses four different branches listed in Table 1. We focus in 
this article on regime 2 which can be associated with the integrable perturbation <p 1 2 of 
the unitary minimal conformal models Ml l+i- The Boltzmann weights in the off-critical 
region are expressed in Jacobi's theta functions depending on a spectral parameter u 
restricted to the interval listed in Table 1, and an elliptic nome < \p\ < 1 which can be 
identified with the ordering field. The latter is of magnetic or thermal type depending 
on whether L is odd or even. As p — > O^ 1 the system approaches criticality. 
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Table 1. The four branches of the dilute A models. Listed are the allowed values of the 
spectral and the crossing parameter entering the Boltzmann weights as well as the central 
charge and the anomalous scaling dimension of the perturbing operator. 

In this article we put forward the leading off-critical excitations in the spectrum 
of the associated transfer matrix in regime 2 for arbitrary L when the thermodynamic 
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limit is taken. This allows us to determine the associated spectral gaps in the bulk limit, 
which determine the finite off-critical correlation length 

r x = - lim ln^H^S^Vbl^r 1 , (p-O*). (1) 

Here A 1; A are eigenvalues of the transfer matrix on a strip of length N correspond- 
ing to the lowest excited and the ground state, respectively. As the critical point is 
approached the spectral gap collapses, and the associated correlation length diverges 
with critical exponent determined by the anomalous scaling dimension A = A 12 of the 
perturbing field 12 . The renormalization group behaviour of the system is encoded 
in the bulk universal scaling functions < S ± which depend on the dimensionless scaling 

1-A' 

variable q'/\p\ 1_A , with q' being a dummy variable for the perturbation by any other 
relevant field with scaling dimension A'. Further information on the renormalization 
group characteristics of the <j) 1 2 perturbation comes from the singular part f s of the free 
energy density of the dilute Al model |3|, f| and its associated bulk scaling functions 
Q ± as p -> ± , 

/ = - lim iV- 1 lnA = / s + / ns , / s ~ \p\ T= ^Q ± (q'/\p\^). (2) 

The explicit form of the bulk scaling functions and Q ± is not accessible via the dilute 
Al model since it realizes perturbation by only one relevant field, i.e. q' = 0. However, 
besides the scaling dimension of the perturbing operator one can extract the following 
combination of critical amplitudes (see e.g. [0]): 

lim/ s .( 2 = Q ± (0) l S ± (0f. (3) 

p^0± 

While each amplitude separately depends on the normalization of the fields and is a 
non-universal quantity, the above combination is universal in the sense that it only 
depends on the universality class determined by the minimal model Ml,l+i an d its 2 
perturbation. Since in the thermodynamic limit the statistical lattice model can be 
effectively described by an integrable quantum field theory the amplitude (|3|) is also 
accessible by exact S-matrix theory. 

The best known example of the exact S-matrices corresponding to the <p 1 2 pertur- 
bation of M l> l+i is Zamolodchikov's construction || related to the simple Lie algebra 
Eg for the Ising model at T = T c in a magnetic field. In the context of the dilute Al 
model this case corresponds to L = 3 in regime 2. Further examples exhibiting the ex- 
ceptional algebraic structures E-j and E 6 are the scattering matrices 0, 0, |5| for L = 4 
and L = 6 which are associated with the leading thermal perturbation of the tricritical 
Ising and tricritical three-state Potts model, respectively (see Table 2). For general L 



it has been proposed by Smirnov [|Kj that the 2 1 and (f) l 2 perturbations should be 
considered as reductions of the Zhiber-Mikhailov-Shabat (ZMS) model [|ll|, |12| . While 
for arbitrary values of the purely imaginary coupling constant the theory violates uni- 
tarity, one can recover physically sensible theories at rational values by reducing the 
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state space. The corresponding scattering amplitudes are obtained at roots of unity, 
q L = 1, from the Izergin-Korepin R-matrix Jl3| which can be associated with the quan- 
tum group U q (A% ). The particle spectrum consists in general of a three-component 
kink and its bound states, called breathers. For L = 3, 4, 6 the aforementioned special 
cases are recovered with L = 3 being distinguished by the absence of any kink state |10| . 
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Table 2. Displayed are three special cases of minimal conformal models together with the 
values of the associated central charge and the anomalous scaling dimension of the perturbing 
operator. All of them are distinguished by an underlying Lie algebraic structure. In case of 
the unrestricted model, L=oo, one recovers the conformal theory associated with a free boson. 

The universal quantity (^]) can be extracted from the exact S-matrix theories by 
applying the thermodynamic Bethe ansatz [113, ITH IIR [T7| which allows the two-particle 



scattering amplitude to be related to the infinite bulk vacuum energy: 

(6) = ^ = -47rm 2 Q ± (0)5 ± (0) 2 , -^mS^) = - f> n e-^l . (4) 
Oi dp * 

n=l 

Here (©) denotes the infinite volume vacuum expectation value of the trace of the en- 
ergy momentum tensor of the associated integrable quantum field theory, m the physical 
mass of the fundamental particle (which can be identified with the spectral gap (|1|)) and 
S(/3) its two-particle scattering amplitude depending on the rapidity (3 parametrizing 
the two-momentum = m(cosh/3, sinh (3). Formula (||) is obtained when computing 
the groundstate energy of the integrable quantum field theory defined through the exact 
S-matrix on a cylinder of circumference R. In the UV limit, R — > 0, the infinite bulk 
vacuum energy can be extracted by comparison against results from conformal perturba- 



tion theory [|Lj, [L6| . Equation (§]) has to be applied with caution as it must be modified 
when logarithmic terms are present in the singular part of the free energy. It is, how- 
ever, believed to hold true for the scaling field theories and S-matrices mentioned and 



the amplitude (f|) has been reported in (T7| based on the exact scattering matrix theory 



in |n|. In addition, from formula @ it can be seen that the universal amplitude enters 



the vacuum expectation value of the energy momentum tensor. The latter is needed as 



important input information in the form factor program [|T^, fL9|| . Setting up a set of 
recursive functional equations it allows all multi-particle form factors of the quantum 
field operators to be constructed starting from their vacuum expectation values. 

Our discussion of the dilute A L models provides an alternative and independent way 
of calculating the universal amplitude P) and is consistent with formula (H) involving 
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the scattering matrices |T0[ and their thermodynamic Bethe ansatz analysis. It is one 
of the main results of our paper. 

We also reveal the close link between the exact S-matrix theory and the transfer 
matrix method on another level which is of practical importance in the calculation of 
the excitation spectrum and the spectral gaps (H). As in the bootstrap approach we 
proceed by constructing functional relations using entirely analytic arguments. In fact, 
we establish a one-to-one match between the functional relations satisfied by the ex- 
citation spectrum of the transfer matrix given in terms of elliptic functions and the 
scattering amplitudes given in terms of trigonometric functions. Besides crossing sym- 
metry and unitarity, the crucial bootstrap identity characterising the S-matrices 
the functional relation associated with a self-fusing process b x b — > b: 



is 



S bb ((3)S bb (P + if) = S bb ((3 + if) . (5) 

Here S bb denotes the two-particle amplitude either of the fundamental particle or of 
a bound state of it. This distinction is important as, except for the case L = 3, the 
fundamental particle does not necessarily display a self-fusing process. However, there 
are always bound states of the fundamental particle which exhibit this peculiar property 
10|| . The equation (|5]) is understood w.r.t. analytic continuation and implies a pole in 



the physical sheet at (5 = if. 

In the context of the dilute Al model we demonstrate in the thermodynamic limit 
that for regime 2 the bootstrap equation (|5|) is reflected in the excitation spectrum of 
the dilute Al transfer matrix by the appearance of an analogous identity, 

r(u)r(u + 2A) = r(u + A) , r(u) — lim t—tt - (6) 

N->oo Aq{U) 

The variable u is the spectral parameter entering the Boltzmann weights and does not 
have the physical interpretation of the rapidity in (^). The constant A is the crossing 
parameter whose dependence on the maximal height value varies in the different branches 
of the model (see Table 1). We stress that the relation (|6|) holds for any excited state A 
not only the leading one. It has been reported before in the literature for L = 3, 4, 6 as 



an outcome of the exact perturbation theory approach ||, |2C| [21], ^ based on various 



string hypotheses |23|, |25], £4], £6| . 

In this paper we derive it in the massive regime for general L directly from the 
known eigenvalue spectrum using only analytic arguments. Functional equations of the 
transfer matrix in the form of an inversion relation are well-known in the literature, 
e.g. 0) |29], |3|. However, we emphasize that the above functional identity (^|) of 



bootstrap type is more powerful as it implies the inversion relation: 

r(u)r(u + 3A) = 1 . (7) 

This relation corresponds to a combination of the quite general properties of unitarity 
and crossing symmetry in exact S-matrix theory, while the appearance of a bootstrap 
equation such as (P) displays a characteristic feature of the model. For L = 3,4, 6, oo 
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we indeed recover all bootstrap identities by employing a generic formulation in terms 



of Coxeter geometry, originally applied in the context of affine Toda field theories [30 



31] , p2| . The fundamental excitations of the dilute ^4^,00 models can be cast into a 
single formula involving the Coxeter element of the underlying Lie algebra, just as the 
ADE affine Toda S-matrices for real coupling have been JJ2], |30|. The latter contain the 
aforementioned scattering matrices as special cases, with the important difference of a 
coupling dependent CDD factor |33|, ||, || [34], |35| . 

As in the context of the bootstrap construction of scattering matrices, the functional 
equations describing the excitation spectrum determine it completely, up to the location 
of its zeroes and poles. The missing information on their position is encoded in the 
Bethe roots. We will argue that in the thermodynamic limit the leading excitations are 
related either to a hole or a two-string in the groundstate distribution of the Bethe roots. 
This picture is supported by revealing the underlying Coxeter geometry related to the 
exceptional algebras £s,7,6 inherent to the special cases L = 3,4,6 and to the simple 
Lie algebra -D4 for the unrestricted model L = 00 (see Table 2). The latter is related 
to the sine-Gordon model at special values of the coupling constant |I(J. Moreover, 



our findings are consistent with all previous numerical investigations in the literature 
23| , |25| , |36| , |37| and the outcome of the exact perturbation theory approach for L = 3, 4, 6 



20, 21, 22 



This article is organized as follows. In Section 2 we briefly review the definition 
and the finite size eigenvalue spectrum of the dilute Al model. We also report a series 
expansion of the free energy density in the elliptic nome based on earlier calculations 
of the groundstate || [| . Section 3 states our analyticity assumptions and presents the 
derivation of the bootstrap equation determining the explicit form of the excitation 
spectrum in the thermodynamic limit up to its poles and zeroes. While the discussion 
in Section 3 applies to any excited state, Section 4 presents a concrete proposal for the 
Bethe root distributions of the leading excitation in regime 2 in the form of a hole and a 
two-string. The position of the zeroes is determined for these excitations and we calcu- 
late the spectral gaps in the eigenspectrum as well as the associated universal amplitude 
ratios. Section 5 compares the outcome against earlier results for L = 3, 4, 6 and presents 
new results on the unrestricted model L = 00. The underlying Lie algebraic structures 
£s,7,6 and -D4 are revealed by using Coxeter geometry which allows the fundamental ex- 
citations of all four models to be expressed in a single formula. We also derive numerous 
new identities which are in one-to-one correspondence with the bootstrap identities of 
affine Toda field theory. The underlying Lie algebraic structures will also be revealed 
in the derivation of the spectral gaps associated with all 8, 7, 6, 4 fundamental particles. 
We obtain a series expansion in terms of the ordering field (the elliptic nome) yielding 
the corrections to the scaling behavior ([!]) in the off-critical regime. The powers in the 
expansion match the affine Lie algebra exponents of £8,7,6 and -D4. Section 6 states our 
conclusions. 
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2 The off-critical dilute Ai model 



The dilute Al model p| is an exactly solvable RSOS model defined on the square 
lattice. Each site is allowed to take one of L possible height values and the associated 
adjacency matrix is given by 3 — A with A being the Cartan matrix of the Lie algebra 
Ai = su(L + 1). In other words, the Boltzmann weights of a face are non-zero only 
if neighbouring heights coincide or differ by ±1. Away from criticality the Boltzmann 
weights are parametrised in terms of Jacobi's theta functions of nome p = e~ T . We will 
frequently make use of their conjugate modulus representation, 



with 

oo 

E(z,p) = ]J(l-p-h)(l-p n /z)(l-p n ) . (9) 

n=l 

The explicit expression of the Boltzmann weights will not matter in our discussion and 
can be found in |J. The eigenspectrum of the row-to-row transfer matrix for a strip of 
length N has been reported in (we choose N even), 

Mu) = »{ WkpWMJ) I G{u) 



+ 



??i(u,p)7?i(3A-u,p) , l JV G(u-X) 



i?i(2A,p)i?i(3A,p) J G(u-2X) 



+ " 1 ^ 1 (2A,^ 1 (3A, P ) / G(U " 3A) (10) 
Here we have introduced for convenience the auxiliary function 

The behaviour of the model depends crucially on the choice of parameters. Henceforth 
we will restrict ourselves to the regime 2 with the interval of the spectral parameter and 
the crossing parameter A as specified in Table 1. The complex numbers u n , n = 1, N 
characterising the eigenvalue are the Bethe roots, subject to the Bethe equations 

^(u+X,p)\ N G(u) 



p(u m ) = -. with p(") = ( fl i(M _ A ; p) ) Giu + x) G(u-xy (12) 

where the phase factor to = exp {i^-A^) with i = 1, . . . , L depends on the particular 
position of the Bethe roots chosen. In fact, using the known transformation properties 
of the theta function 

(u ± ir, e~ r ) = -e ± ^e^ 2iu -& 1 (u,e- T ) , r>0 (13) 

$1 (u + 7T,p) = — l & 1 (u,p)='&i(—U,p) (14) 
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one verifies that the corresponding eigenvalue A(u) is invariant under the simultaneous 
replacement 

u n — ► u n ± it and uj — ► we T ' 4 (15) 

as well as under the shift 

U n ^U n ±7T (16) 

of a single Bethe root. Moreover, one easily deduces from the transformation properties 

G(u) =e ±miX G(u±tr) =G(u + n) (17) 

of the auxiliary function ( |TT] ) that the eigenvalues satisfy 

A(u + tt) = A(u) and A(m ± ir) = ]9 TAf e ±i6AAf e T4A ^A(w) (18) 

while the function determining the Bethe equations (|12|) is doubly periodic 

p(u) = p(u + ir) = p(u + 7r) . (19) 

Finally, from the identity i?i (w, e 4?r p) = e l 5i?i it is seen that, for all finite N, 

the eigenspectrum and Bethe equations are invariant under the exchange p — > e l7T p as 
only ratios of Jacobi's theta functions occur in (|T0"D and (i~2|). However, the Boltzmann 
weights H are in general not invariant under sign reversal of the elliptic nome, but 
can be related through a relabelling of the heights if L is odd. Therefore, while all 
eigenvalues of the transfer matrix can be parametrized through fllQl) and (|12|) with p 
replaced by \p\, an eigenvalue present for p > might be absent for p < 0, and vice 
versa, if L is even. 



2.1 The groundstate and scaling corrections 

The eigenvalue A corresponding to the groundstate of the model has been investigated 
for arbitrary L in [|J by applying Baxter's exact perturbation theory approach |38| in the 
ordered limit \p\ — > 1. The groundstate Bethe roots are assumed to lie on the imaginary 
axis, which has been supported by numerical investigations. The result obtained for the 
groundstate is assumed to hold for any < \p\ < 1 in the thermodynamic limit and 
allows one to compute the free energy density (0) 



/ = - lim A- 1 lnA (f ) = 4^ V 



COS 



'5tt 



-g)(6^+l) 



V 



&(6fc+l) 



£ (6fc + 1) sin (6fc + 1)] l-ps 



-(6fc+l) 



cos 



'5tt 



-(6fc-l) 



^ (6A; - 1) sin [g (6A: - 1)] 1 - p^ (6fc - 1} 



oo 



k=l 



sin 2 ^fccos5A/ccosA/c p 



2/,- 



k 



cos 3AA; 



P 



2A- 



fs fns- 



(20) 
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This formula is obtained from the result for the partition function per site derived 
in || |J by applying the Poisson summation formula and the residue theorem. The 
spectral parameter has been set to the isotropic value u = 3A/2 in order to ensure that 
the eigenvalues are real. We come back to this point later on. From the series expansion 
(p0|) the coefficient of the leading term in the critical limit p — > + is obtained as 

, , s l r- cos t?t(5A — 7r) , x 

Q+(0) = lim p~—f s = 4V3 ^ 21 

p^o+ sin ^7r 

with 2 A = 2Ai 2 = 2 — 6X/tt being the scaling dimension of the perturbing operator. As 
already pointed out in the introduction, the latter depends on the chosen normalization 
of the fields and is a non-universal quantity. In order to obtain the universal quantity 
(|D we need to compute ([D from the leading excitation in the thermodynamic limit. To 
this end, in the next section we derive functional equations satisfied by any excited state 
in the thermodynamic limit. 



3 Bootstrap identities in the dilute Ai model 

Conventionally the excitation spectrum in the thermodynamic limit is calculated pro- 
ceeding via the thermodynamic Bethe ansatz. Alternatively, one might proceed entirely 
by analytic arguments as has been demonstrated e.g. in the context of the XYZ model 
|2"9fl . Rather than solving the intricate non-linear integral equations of the thermody- 



namic Bethe ansatz, we will instead derive a set of functional relations for the ratio 

r(u) = lim — Vv • 22 

V ' N^oo Aq(u) 

Here A(u) denotes the eigenvalue of any excited state over the ground state A (u) such 
that the ratio A/A has a non- vanishing well-defined limit as N — > oo. We will in par- 
ticular show the functional equation (^), which, in conjunction with certain analyticity 
assumptions, imposes powerful restrictions determining the excitation spectrum up to 
the position of its poles and zeroes. The latter are fixed by the Bethe roots which 
are usually unknown. However, the functional relations we are going to discuss yield 
constraints on their positions. 



3.1 Asymptotic behaviour in the thermodynamic limit 

Recall that the Bethe roots belonging to the groundstate Ao lie on the imaginary 

axis. It is reasonable to assume in general that for the low lying excited states A the 
associated Bethe root distribution {u n } is similar, and that at most a finite number of 
holes are present or a finite number of Bethe roots have non- vanishing real part. Well 
away from the critical point, < | hxp\ <C 1, this information is sufficient to show that as 
N — > oo either the first, second or third term in the eigenvalue expression ( |10|) becomes 
dominant as the spectral parameter sweeps through the allowed interval < u < 3A. 
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The crucial feature exploited is the conjugate modulus transformation (||), as outlined 
in Appendix A. Explicitly, one finds for the groundstate as well as for the excited states 
the following: 



A(u) ~ < 



i9i(2A-m))?i(3A--u) 
th(2A)0i(3A) 



N 



o(e 



-N 



G{u) [I 

i?i(u)tfi(3A-u) T" G(u-A) N , n („-N\l 
i?i(2A)i?i(3A) f G(u-2A) L 1 ^ U v C /J 

G(u-3A) + 



(2A)i?i(3A) 

I I ^i(tt)i?i(A-u) 
w \ i?i(2A))?i(3A) 



A' 



o e 



< u < X 
X < u < 2A 
2A < u < 3A 



(23) 



By the above notation we indicate that all terms besides the leading one are exponen- 
tially suppressed in N. In order to discuss the thermodynamic limit we introduce 



9{u) 



lim g N (u) 



9 N (u) 



G(u) 



(24) 



where Gq(u) denotes the auxiliary function ( p|) corresponding to the groundstate. By 
an analogous argument as in the XYZ case ||29|| , the convergence of the above limit is 
expected to be uniform as only a finite number of Bethe roots have non-vanishing real 
part or are missing compared to the groundstate distribution. The excitation spectrum 
then obeys 



rlu) 



lim 



A(«) 



9W 



< Re u < A 



g(u - X)/g{u - 2A) , A < Reu < 2A 
^giu-SX)- 1 , 2X<Reu<3X 



N^oo A (u) 

An obvious consequence of this limiting behaviour is the identity 

r{u)r{u + 2A) = r(u + A) , < Reu < X. 



(25) 



(26) 



In order to promote this relation to the functional equation (|B]) valid in the whole 
complex plane, we need to investigate the analytic behaviour of the excitation spectrum 
in the thermodynamic limit. 



3.2 Analyticity and the general form of the excitation spec- 
trum 

From the explicit expression flT0p for the eigenspectrum, it is inferred that, for finite N, 
poles and zeroes of the ratio A(u)/Ao(m) are completely determined by the location of 
the Bethe roots. Employing the asymptotic behaviour (^) one easily deduces that for 
N ^> 1 the purely imaginary Bethe roots u n of the excited and of the ground state give 
rise to poles located at 

uf m = u n + A, u n + 2A mod7r . (27) 

These pole contributions will cancel each other when the limit N — > oo is taken, since 
the imaginary roots then become densely located along the imaginary axis centered at 
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A, 2Amod7r. The only exceptions are poles in the groundstate which give rise to zeroes 
of the eigenvalue ratio A(m)/A (m) when holes are present in the excited state. Bethe 
roots u n of the excited state belonging to a string, i.e. Re-u„ 7^ 0, will also give rise to 
isolated zeroes and poles. However, since the number of holes and strings is assumed to 
be finite, one infers from (Pp that in each of the intervals 

0<Re?i<A, A<Rew<2A and 2A<Reu<3A (28) 

the ratio A(m)/A (u) is analytic except for isolated points, and therefore converges to- 
wards a meromorphic function r(u) = liniTv^oo A(u)/A (w) in these intervals. Upon 
analytic continuation r(u) can then be extended to the whole complex plane and the 
identity fl26|) becomes the functional equation (|6|) valid for all complex values of the spec- 
tral parameter. As an immediate result we have the inversion relation (]?]). Applying 
the inversion relation twice we deduce that the excitation spectrum is doubly periodic: 

r(u) = r(u + it) = r(u + 6A) . (29) 

The first period is a direct consequence of (|l"8"D and is already present at finite N. The 
assumption of analyticity together with the double-periodicity property are powerful 
restrictions which determine the excitation spectrum up to the position of its poles and 
zeroes by applying the following theorem |38|| : 



Theorem. Let f(u) be a meromorphic function satisfying 

f(u + 2K) = f(u + 2iK') = f(u) 

with K,K' being quarter periods to the elliptic modulus k. If f(u) has n poles per 
period rectangle then it has also n zeroes and can be expressed as 









-uf) 




^'1 


01 




u 


-0 


,e~ 


„ K ' ' 



f(u) = Ce^H 



Here C is a constant and s' an integer determined by the sum rule 

n 

J2( u f - u f ) = 2sK - 2is>K ' ' s,s'eZ. (30) 

3=1 



We therefore conclude that the excitation spectrum of the dilute model at infinite 
volume is of the form 
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Inserting this expression in the functional relations (0) and (||) we deduce several con- 
straints on the location of poles and zeroes. From the inversion relation it follows that 

C 2 = 1, s ' = 0, s = 3An and uf = uf ± 3A . (32) 

It is therefore sufficient to locate either the zeroes or the poles. From the functional 
equation (|6|) we deduce that the existence of a zero uf implies another zero either at 
uf + A or at uf — A. Thus, all zeroes come in pairs and n in ([!!]) is even. Furthermore, 
the constant in ([H]) is fixed to be C = (— l)"/ 2 . 



3.3 Hermitian analyticity and unitarity 

In order to restrict the position of zeroes and poles further we derive another functional 
equation which is the analog of hermitian analyticity and unitarity in exact S-matrix 



theory [39 



r(-u*Y = r(w)- 1 . (33) 



Here * denotes complex conjugation. In order to prove Q33|) we need to make the further 
assumption that the complex Bethe roots form a set which is symmetric w.r.t. to the 
imaginary axis modulo multiples of tt, i.e. the set {«„} should be invariant under the 
mapping 

u n — > — u* n mod7r . (34) 

This is trivially fullfilled for the Bethe roots of the groundstate Ao, which lie on the 
imaginary axis. In the case of an excited state A a finite number of the Bethe roots 
{u n } might have non-vanishing real part and the above constraint is non-trivial. However 
there is numerical evidence for L = 3, 4 that ([54]) is satisfied and all string hypotheses 



23] , Pq , pq , p2| previously reported in the literature for L = 3, 4 are subject to this 
constraint. As a consequence of ( |3"4"P the auxiliary function ( |TT]) for the groundstate as 
well as for an excited state satisfies 

G(-uT = TT + = GiuT 1 ■ (35) 

From this property we infer for all N that the corresponding eigenvalues QH^ ) are crossing 
symmetric: 

A{u) = A(3A-m*)* . (36) 

Upon analytic continuation in the thermodynamic limit we conclude from this crossing 
symmetry and the inversion relation (0) that the property (|33"D must hold. Thus, we 
can eliminate the pole variables u p from (|31[ ) by setting u p = —u z *. We further deduce 
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from crossing symmetry that the set of zeroes of the excitation spectrum fl3~T|) must be 
invariant under the mapping 



u z ^3X-u z *. (37) 

Together with our earlier observation that the zeroes come at least in pairs (u z , u z + A) 
this now implies that one has either a pair 



{u z ,u z + X) with Reu z = X (38) 



or a quartet of zeroes 



u z ,u z + X,2X-u z *,3X-u z *) . (39) 



As argued earlier, the real parts of the zeroes u z are determined by the complex Bethe 
roots u n with Rew„ ^ 0, whence the conditions (|38| ) and ( |39"D constitute further restric- 
tions on the string structure of the dilute Al model, besides (|34j). We leave a complete 
investigation of the string structure to future work f4(J and now focus on the leading 
excitations in regime 2. 



4 The leading excitations in regime 2 

Having determined the general form of the excitation spectrum in regime 2 in the ther- 
modynamic limit, we now proceed by putting forward a concrete proposal for the leading 
excitations which will ultimately allow us to compute the leading correlation length ([!]) 
and the universal quantity (H). 

Proposal. The leading excitation in regime 2 consists either of a hole or the presence 
of a two-string in the groundstate distribution of the Bethe roots. 

This proposal is supported by numerical calculations [^3|, |25|, |36], [J7|, the thermo- 
dynamic Bethe ansatz p3f and the exact perturbation theory approach for the dilute 



Al models in regime 2 with L = 3, 4 |, |20], [H], p2"| . Further support comes from ex- 
tracting the correct Lie algebraic structures for L = 3, 4, 6, oo, as we will discuss in the 
subsequent section. Here we first discuss the implications of our proposal for general L. 

From ( p3[) and fl2~5|) we infer that the zeroes of the excitation spectrum are determined 
by the zeroes of the auxiliary function (p4|) in the appropriate domain of analyticity. 
The latter are determined by the Bethe roots according to the definition (|TT|). In the 
thermodynamic limit only the real part of the zeroes u z in the excitation spectrum (|3~TD 
needs to be determined, since lmu z is linked to the imaginary part of the Bethe roots 
which can be parametrized by a continuous variable as iV —>■ oo. Given a Bethe root we 
note that G(u) in general has a simple zero and pole at 

G(u) : uf = u n — A and uf = u n + A, (40) 
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respectively. Accordingly, the second and third terms in the eigenvalue expression fllCf ) 
exhibit zeroes and poles located at 



G{u - A) 



ufj = u n , u n + 3A , ufj = u n + A, u n + 2A (41) 



G(u - 2A) ' 
and 

Giu-SX)- 1 : uf II = u n + AX, u? n = u n + 2\. (42) 

We deduce from the asymptotic behaviour (|23|) that in the thermodynamic limit these 
poles and zeroes are only of significance for the excitation spectrum if their real parts 
are located inside the domains of analyticity stated in (|25|). This constraint originates in 
the fact that the excitation spectrum ([31]) is the analytic continuation of meromorphic 
functions whose region of analyticity is restricted. The resulting zeroes are therefore 
characteristic for the regime chosen. Recall from (|T7|) that the zeroes and poles of the 
auxiliary function are only fixed modulo tt. The following relation for regime 2 will 
therefore prove important to determine the zeroes and poles of the excitation spectrum: 

regime 2 : n = 3A + - , r = A7i = 4^^. (43) 

r L — 2 

Henceforth, it will always be understood that zeroes related by multiples of 7r can be 
identified. 

4.0.1 Zeroes from imaginary Bethe roots and holes 

Before we turn to complex Bethe roots with nonvanishing real part, we first investigate 
the possible zeroes and poles related to the imaginary roots in the thermodynamic limit. 
As we have established in (|40D, ( |41~D and fl42|), only the poles u n + A, u n + 2A may occur 
as N — > oo. Note that these poles are located at the border of the domains of analyticity 
but are common to both relevant terms in (^). These poles cancel against the ones from 
the groundstate function unless there is a hole in the excited state in the thermodynamic 
limit. Then we obtain zeroes with real part located at 

Reu z = X,2X (44) 



in accordance with our earlier observation fl38l) . Whether or not a hole is present can 



in principle be inferred from the Bethe equations, which we will discuss below in the 
ordered limit. In general a hole might be accompanied by a string, i.e. some of the Bethe 
roots might in addition have vanishing real part. This picture has been suggested by 



unpublished numerical calculations [26] for L = 4 where the presence of a hole excitation 



has been reported in connection with the following 3-string: 

v,v±2X, Rew = f. (45) 



This scenario has been investigated by the exact perturbation theory approach ]22] and 
it was observed that the three string does not of itself contain any relevant information. 
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We confirm this outcome in our analytic approach; from (f40|), (f4T|), and ( |4*2| ) we deduce 
for general L that the proposed 3-string does not give rise to any zeroes in the excitation 
spectrum. 

Further restrictions on the existence of a hole excitation might originate in the max- 
imal height value L and in the choice p > or p < 0. In fact, when we evaluate these 
zeroes for L = 3, 4, 6, oo in a subsequent section we will see that it does not give the 
leading excitation for L = 3. This is consistent with the suggested picture from scat- 
tering theory [10] if we associate the hole-excitation with a kink state, absent for the 
Ising model in a magnetic field. The fundamental particle in the case when the hole 
excitation is not present should be associated with a particular two-string, which we 
consider next. 



4.0.2 Zeroes induced by a two-string 

A general two-string subject to the constraint (|3"4] ) is of the form 

v ± x, Re v = or ? , Im x = . (46) 

Setting x = A/2 and Rev = it/ 2 this type of string has been observed numerically for 
L = 3,4 |^3[ ^6|. We propose it here for general L and show below that it leads to 
the correct algebraic structure for the special cases L = 3,4, 6, oo. In comparison with 
the exact S-matrix theory proposed in [IIJ , the associated excited state of the dilute 
model should be identified with the lowest breather, i.e. a bound state of two kinks. 
One easily verifies that G(u) has simple zeroes and poles at 

uf = v-f,v-± and u? = v + %,v + ^. (47) 

Since Rev = n/2 the second zero has real part greater than A and thus does not 
constitute a zero of the eigenvalue in the thermodynamic limit. Note also that the poles 
are not present in the thermodynamic limit. The second term and third terms in the 
eigenvalue expression fllQ~D exhibit zeroes and poles at 

ufj = v-§,v + ±,v-l-^ v + l-z (48) 
uf t = v + l,(v + f) 2 ,v-l-^ 

and 

ufn = v + ±-Z,v + f and u f n = v + f ,v + f . (49) 

We therefore deduce the following location of the real parts of the zeroes in the excitation 
spectrum when the thermodynamic limit is taken: 

Reu z = |:,A + £,2A-£,3A-£. (50) 

Note that the resulting quartet fulfills the condition (|3~9p. 

Having determined the zeroes of the excitation spectrum we are now in the position 
to write down the corresponding ratio of eigenvalues ( |31] ) and to calculate the resulting 
spectral gaps. Before we come to this point we further support our proposal on the 
excitations and their Bethe root structure by investigating the Bethe equations (|T2|) in 
the ordered limit. 
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4.1 Bethe equations in the ordered limit 

As solving the Bethe equations directly is a complicated matter, we consider them in 
the ordered limit \p\ — ► 1 where we can check the proposed excitations for consistency 
by investigating the asymptotic behavior as iV — > oo. We will show that the function 
(H^) characterizing the Bethe equations is a pure phase in the massive regime for N ^> 1 
as required. Its asymptotic behavior is determined by the real part of the Bethe roots 
which is the essential information entering the excitation spectrum (|3ip . The latter 
depends continuously on the elliptic nome, whence we expect that our findings in the 
ordered limit give a qualitative picture valid for all values of \p\. 

As we did for the eigenspectrum of the transfer matrix in Appendix A, it is also of ad- 
vantage here to rewrite the Bethe equations using the conjugate modulus transformation 



^ j E(x 2s w) \ N yr g E(x 2s w/w n )E(x 4s w n /w) 
P[W) ~ \ wE(x 2s /w) J l} x Wn E(x 2s w n /w)E(x is w/w n ) ' 1 j 

We recall the definition of the variables, w = e~~?~ u , s = ^| and x = e - ^ with 
r = . The dependence on the elliptic nome x 2r is suppressed in the notation. 

4.1.1 Bethe equations and the occurrence of holes 

Let us start by investigating the ordered limit x — > (\p\ — ► 1) when only purely 
imaginary Bethe roots are present. From the conjugate modulus representation (|5T| ) we 
infer that the leading contribution for any Bethe root w m = e ~ Um is given by 

N 

P(w m ) ~ (-Wm)-" l[wf =0(1) . (52) 
n=l 

Thus, p(w m ) lies on the unit circle as expected. Since we are only interested in the 
thermodynamic limit, the exact position of the Bethe roots does not matter. But we 
can extract an additional piece of qualitative information, namely that the N Bethe 
roots w m are completely determined by an equation of order N. This shows that holes 
are absent for any finite iV when only imaginary roots are present. The situation changes 
when we consider in addition Bethe roots which have no n- vanishing real part, since they 
alter the asymptotic behavior in the ordered limit. 

We have mentioned before the possible occurrence of the 3-string (|I5|) for L = 4 
|2T| P% which according to our analysis does not give rise to any zeroes in the excitation 



spectrum for general L. Its significance is the induction of a hole in the purely imaginary 



Bethe root distribution as observed in |22fl for L = 4. We briefly review the argument 
for general L. Set = zx 3s+1 , u>at_i = z'x~ s+1 , = z"x 7s+1 with z,z',z" being 

pure phases up to exponentially small corrections in N. As the remaining Bethe roots 
w n ,n < N — 2 also lie on the unit circle, one deduces from ( |5T| ) the following leading 
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term in the ordered limit (N even ^> 1), 

(zz'z")^ 1 ^ 2* 

p(w m ) ~ w- N+2 ± ^— \{ w n r = o(l) , 1 < m < N - 2 (53) 

^ n=l 

As required p(w m ) lies on the unit circle. In contrast to ([52]), however, the N — 3 
imaginary Bethe roots are now determined by an equation of order N — 2. The extra 
solution points towards the existence of a hole. In addition to the Bethe equations (|53|) 
one also needs to verify that the equations for the complex Bethe roots wn, wn-2 
are of order one in the ordered limit. The product of the Bethe equations corresponding 
to the three roots in the string has to be taken in order to avoid complications in 
the thermodynamic limit where terms either in the numerator or denominator become 
small. Again one finds that the leading term lies on the unit circle in the ordered limit, 
p(w N-2)p{w n-i)P( w n) ~ °(1)> showing that the string structure is consistent. 



4.1.2 Bethe equations for the two-string excitation 

We now consider the Bethe roots for the excitation associated with the two-string. We 
choose the first N — 2 Bethe roots w n to lie on the unit circle and set wn-i = zx 2s+1 , 
wn = z'x 4s+1 . In the thermodynamic limit we have \z\ = \z'\ = 1 and z = z', thus for 
sufficiently large N we may assume these properties up to exponentially small corrections 
in N. For the Bethe roots w m , m < N — 1 we find (6s + 2 = 2r) 



7V-2 



2 s 



p(w m ) ~ {-w rn f- N z^- 2 n Wn = oil) , X 0, N > 1 (54) 

n=l 

while the leading term for the case of the complex Bethe roots reads 

Af-2 

piwN-JptwN) ~ (zz')^- 2 ]Jwf~ 2 = o(l) , x^0,N^>l. (55) 



n=l 



Thus we find for the two-string also that the arrangement of the real parts of the Bethe 
roots is consistent with the Bethe equations in the ordered limit. Note the absence of 
any holes. As indicated before, this is a preliminary check on the string structure and a 
more detailed analysis will be presented elsewhere EH . 



4.2 Spectral gaps of the leading excitations 

We now insert our results (f44|) and ([50]) on the zeroes of the excitation spectrum into the 
general formula (|3i~|) . We recall that the imaginary part of the zeroes is undetermined, 
but in the thermodynamic limit becomes a continuous variable, say lmu z = v G R, 
which parametrizes a band of eigenvalues corresponding to the excited state. In the 
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field theoretic context it can be identified with the rapidity. For the hole-excitation we 
then find 



r ( q .\ Mrx-^-i\p\^)Mrx'^-h\p\^) r , fil 

In general the eigenvalues in the excitation band will be complex and in order to deter- 
mine the spectral gap, or equivalently the correlation length, one has to integrate over 
the whole band. However, from the crossing relation (|36|) we infer that the eigenvalues 
are real at the isotropic point u = 3 A/2. In addition, we have to tune the parameter v 
such that we reach the bottom of the band. The latter value is derived to be v = 
and the correlation length is then obtained via the formula 

4hoie(P) - -lnr hole ( T ) 1 , =f& ~ 2[n ^^ E _ E • (■■>') 



6A 



Using a standard identity for the elliptic functions we obtain the series expansion in the 
elliptic nome 



= 4v ^ ^Ti i - |p|&(^) + 473 ^ 6^1 i - |p|^(6»-i) ' ^ 

n=0 ,Jr 1 n=l 1 



from which we easily extract the amplitude (H) of the leading term in the critical limit 
p - 0± 

^(O)- 1 = lim IpI-^^Cp)" 1 = 4v/3 . (59) 



Together with the general expression for the singular part of the free energy ( pT|) this 
puts us in a position to derive the universal amplitude ratio (H): 



sin? L 



C ± (0)5 ± (0) a = JT^fW , ^>3- (60) 

This coincides with the result from the thermodynamic Bethe ansatz analysis re- 
ported in [|17]] based on the exact scattering matrices constructed in [10 when the fun- 



damental particle is related to a kink. It needs to be pointed out that the general result 
agrees with the outcome for the special cases L = 4, 6 previously reported in p0| . 



For L = 4, i.e. for the leading thermal perturbation of the tricritical Ising model, (|60| ) 
can also be compared against the numerical findings obtained in JII| . 



The case L = 3 is special as kinks are absent from the particle spectrum |IJ. In 
the context of the dilute A 3 model, also, numerical investigations |23|, of the Bethe 
equations did not show evidence for the presence of a hole. As we have pointed out 
before, further exceptions to the hole constituting the leading excitation might occur 
for L even, when the sign of the elliptic nome is negated. This is, for example, the case 
when L = 4. 
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For these exceptional cases we therefore consider as the leading excitation the two- 
string originally proposed in |23] for L = 3. In the context of the scattering theory flO 



this type of excitation would correspond to the lowest breather, which is the lightest 



particle exhibiting a self-fusing process. From ( pTD and (|50|) we calculate the associated 



spectral gap by repeating the same steps as before and find 



SstringlPJ - „Q fn nL-2 |„|^.\„q ( n , tt 3L+2 i f 



E 



*V4 12L+2'I^I y^4V 4 i 12 L+2 ' ^ I 

i ^ — — (sin ^¥r# + sin 

i n i-bl 5 " 

7T 3L + 2 , • 7T 

6 L+2 ' 0111 6 L+2 



L+2 1 ° 6 L+2- 



(sinfff + sin£fef)bK |phO, (61) 



We point out that the sum (^H]) runs at most over the integers relatively prime to 6, i.e. 
n — 1, 5, 7, 11, ... These integers coincide with the allowed spins of integrals of motion 
compatible with a self-fusing process |7| . In order to obtain the universal amplitude (Q) 
in the case that the hole excitation is absent we need only to determine the ratio of the 
two different correlation lengths in the critical limit, 



£string(0) _ Sm 3 



WO) cos + sin |y 



< 1 , L > 3. (62) 



As expected, we deduce that in presence of both a hole and a two-string the smallest 
spectral gap is given by the hole-excitation. In the context of the scattering matrices 
constructed in [|l(J expression (|6^) determines the mass ratio of the kink and the lowest 
breather. 



5 The excitation spectrum for L = 3, 4, 6, oo 



For L = 3, 4, 6 and regime 2 the excitation spectrum has been previously put forward in 
the literature case-by-case, using the exact perturbation theory approach in the ordered 
limit 0, ^(J, |22| and the different string hypotheses based on numerical data p3[ [25|, 
i, |26|, [36j . A significant role in these proposals is played by the underlying Lie algebraic 
structures related to the exceptional algebras E 8 j^ which also arise in the associated 
exact scattering theories. 

The excitation spectrum for the unrestricted model (L = oo) connected to the alge- 
bra D4 has not been previously considered and its formulation is a new result. 

We show that our proposal of a hole-excitation and a two-string excitation leads 
to the correct Lie algebraic structures by connecting the location of the real part of 
the zeroes (fB|) and (50) to the Coxeter geometry of the respective algebras. Invoking 
similar Lie algebraic Coxeter identities as were originally applied in the context of affine 
Toda field theory [|3(], [32]] , we summarize the previously reported excitation spectra for 
L = 3, 4, 6 in a single formula and prove additional functional equations besides the 
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ones already stated for general L. These relations completely fix the location of the 



zeroes in (]31f) for the higher fundamental excitations J23|, |36], |, |20], [H], as well, 
and can be regarded as equivalent to a string hypothesis. We start by stating the 
excitation spectrum in generic Lie algebraic terms for L = 3, 4, 6, oo and then identify 
the excitations belonging to the hole and the two-string. 



5.1 Excitations in terms of Coxeter geometry 

In the following let q stand for one of the simple Lie algebras E§j^ or D^. Then 
we expect, depending on the sign of the elliptic nome, at most rankg fundamental 
excitations or eigenvalue bands in the thermodynamic limit. Furthermore, we denote by 
h the Coxeter number of q which obeys the following identity in terms of the maximal 
height value, 

k=6§±§. (63) 

Introducing the simple roots and the fundamental co-weights A^ of the Lie algebra g 
we define a matrix of functions \i : Z — > |Z by setting 

/^■(2a - 2ft) = -| (AJ, a a Oi) , c t = ±1 . (64) 

Here a denotes the Coxeter element of the Weyl group of q. It is uniquely defined 
w.r.t. to a specific bicolouration of the Dynkin diagram specified by the "colours" Cj. 
See Appendix B for details. We are now prepared to state the excitation spectrum for 
L = 3,4, 6, oo: 

r Au) = (-l) A *i If < — 7 f > • 65 



0=1 



Here j = 1,2, ...rankg runs over the fundamental particle excitations and the index * 
specifies the root in the Dynkin diagram of the simple Lie algebra g which is connected 
to the affine root in the Dynkin diagram of the affine extension [2~I]. We will 



show below that this particular node is singled out by the two-string. The variable 
a G M parametrizes the imaginary part of the Bethe roots in the thermodynamic limit 
and scales the eigenvalue bands. The sign factor in front is determined by the inverse 
Cartan matrix which is given in terms of the simple roots as = 2 (a^, aci) / (afc, 
In fact, the Cartan matrix allows for an alternative definition of the exponent function 
(|64T), which will be of great practical use in the subsequent computations. Consider the 
g-deformed Cartan matrix [A] q defined by 

(M,)«:=[Mn [< := 7^T- ( 66 ) 
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Then the functions ( |64"D can be implicitly defined by the relation [|^, [13 



i 2h ^/l 

A a=l 

The above relation is non-trivial as it states that the inverse g-deformed Cartan matrix, 
which generically is a rational function in q, simplifies to a polynomial when multiplied 
with the factor (1 — q 2h ). The functions /i then simply count how often the monomial 
q a occurs in this expansion. That the two definitions of \i are indeed equivalent has 
been proven in |43|| . One can also now check explicitly that the above expression 
(|B"5]) correctly reproduces the case-by-case results derived previously ||, ^UL EH], |2"2]1|*|. 



There and in ]36| it was noticed that the characteristic integers a fixing the real part of 
the zeroes in the excitation spectrum flGUP match the integers appearing in the building 
blocks of hyperbolic functions used in describing the ADE afline Toda S-matrices ||. In 
the context of affine Toda field theory the relation of these integers to Coxeter geometry 
was established in |32| . 

One can now explicitly check that the resulting zeroes (|50f) and fl44]) from the hole and 
two-string excitation match the Lie algebraic formula when the node * is specified 
by the simple root connected to the affine root of g. Employing the helpful identities 

A = — - and 2r = 2Ari = h + 2 
2r 3 1,2 

one finds the real parts of the zeroes listed in Table 3 (a) and 3 (b). We note that the 
hole gives always the lightest particle in the mass spectrum (except for L = 3) while the 
two-string associated with the node * corresponds to the lightest particle which exhibits 
a self- fusing process (|5]). It is remarkable that from the simple excitation of a hole and 
a two-string one already sees the four different Lie algebraic structures emerging. 



L 


algebra 


4±2 Reu z 


4 


E 7 


6,12 


6 


E 6 


4,8 


oo 




2,4 



(a) 



L 


algebra 


-k 


^ Reu z 

TT 


3 


£ 8 


1 


1,11,19,29 


4 


E 7 


2 


1,7,11,17 


6 


Eq 


2 


1,5,7,11 


oo 


D± 


2 


1,3 2 ,5 



(b) 



Table 3 (a) and 3 (b). Real parts of the zeroes in units describing the excitation 
spectrum corresponding to the presence of a hole (a) and the two-string (b). The results are in 
accordance with the general formula (^) and display the Coxeter geometry of the underlying 
Lie algebra. Note that the hole excitation corresponding to the kink is absent for L=3. For 
L>3 this excitation corresponds to the particle index j=l. The numeration of the nodes is in 
accordance with the conventions in ||, |2^, 21, 

* Setting obtain the leading contribution previously reported in the literature. 
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The excitations fl65|) with particle indices j > 2 ought to belong to heavier particles 
which can be expressed as bound states of the former ones. This suggest that the 
corresponding Bethe roots distributions, which might contain additional strings or holes, 
have to fulfill further constraints originating from bootstrap identities other than @. 
At the moment we leave further exploration of the string structure of the dilute Al 
models to future work |4(| and simply state these conjectured excitations as additional 
consistent solutions to the equation @ which naturally arise from the Lie algebraic 
structure. Nevertheless, as a preparatory step for a closer investigation of the string 
structure, we state the additional bootstrap identities arising from Coxeter geometry in 
the next section. 



5.1.1 Bootstrap identities 

We start by verifying that the explicit expression ( j65| ) for L = 3, 4, 6, oo obeys the func- 
tional relations (0), (|36D and ([|). To this end it is helpful to introduce the meromorphic 
function 

t(u,a):= ) f (68) 

+ if. w*) 

which enjoys the following transformation properties 

t(u+6X, a) = t(u,a ± 2h) = t(— m, a) -1 = t(u, —a) -1 = t(u,a) (69) 
t(u + 3X,a) = -t{u } h-a)~ l (70) 
t(u + ir,a) = e 2m U(u,a) . (71) 

Using the above identities the verification of the functional equations (|33|), (|36"D and (|6]) 
becomes a purely algebraic problem which can be formulated entirely in terms of the 
generating function (B4I). The latter has been thoroughly investigated in the context 



of affine Toda field theory [^, and proceeding in an analogous manner one can 

establish a one-to-one match between the functional relations satisfied by the excita- 
tion spectrum of the dilute ^3,4,6,00 model and the functional identities obeyed by the 
factorizable scattering matrices of the corresponding affine Toda field theory. 

The properties of the generating function (|64"D which we exploit are p2|, [42], [43| 



/i(a) = n(2h + a) = -fi(2h - a) . (72) 



An additional remarkable identity which is not straightforward to prove reads p2l H 



^ki{h- a ) = Wa( a \ ( 73 ) 

where k denotes the node in the Dynkin diagram of the Lie algebra g which is obtained by 
applying a possible Dynkin diagram automorphism to the node k. The second definition 
of the generating function ( p7| ) is most useful for verifying the formulas 

/i(a) = n(a)* and //(a + 1) + n(a - 1) = (2 - A)fi(a) . (74) 
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We note that the relation ( |33[ ) corresponding to hermitian analyticity in exact S- 
matrix theory is trivially satisfied due to (|S9|). The inversion (|7D and crossing relation 
(|36|) then follow by observing that the preferred node * is self-conjugate, i.e. * — and 
employing (|70|) . The number of zeroes is given by the inverse Cartan matrix element 
A^, which in the case of the E-type algebras coincides with the Kac or Dynkin labels: 

a=l a=l 

From this algebraic identity together with ( |7T| ) we now also read off the double-periodicity 



(p^) of the excitation spectrum (|65|). 

The central equation (|5|) requires a deeper algebraic analysis using Coxeter geometry. 
In fact, it is a special case of more general identities, which, along the lines of the 
analogous analysis in affine Toda field theory |$0|, [32(1 , can be formulated as follows (see 
Appendix B on Coxeter geometry for the explanation of the various symbols): 

Whenever there exist three elements in the Coxeter orbits Qi, associated with the 
simple roots j t = ciat, c\ = ±1 and I = i, j, k which sum to zero, i.e. 

<tH = & Yl V(a ± 7j,) = , Vl := -2£, 

l=i,j,k l=i,j,k 

for some triplet of integers one has the identity 

J] r { (u + 3A|) = l. (76) 

l=i,j,k 

For i = j = k we recover equation (|G]) setting = 0, 2h/3, 4h/3. In the context of affine 
Toda field theory these "bootstrap" identities occur when an intermediate bound state 
is formed, i + j — ► k with k denoting the antiparticle of k. We see that in general these 
identities involve different particle indices. A few examples are listed in Table 4. 



L 


Coxeter identity 


bootstrap equation 


3 


a 2 + cr 13 «3 + a 23 ax = 


r 2 (u + f A)r 3 (n + 2A)r x (u) = 1 


4 


a 5 + a 7 a 7 + a l4 ai = 


r 5(^ + t X)r 7 (u + |A)r 4 (n) = 1 


6 


«i + cr 5 «3 + a w a 2 = 


ri(« + 5A)r 3 (w + |A)r 2 («) = 1 



Table 4. Displayed are three bootstrap identities for the excitations of the dilute ^3,4,6 models 
and the corresponding Coxeter identities in the root systems of the associated £?8,7,6 algebras. 
The numeration of the nodes is in accordance with the convention in [Q, 20, 22]. 

In the present context, these bootstrap identities completely fix the position of the 
poles and zeroes and are therefore equivalent to a string hypothesis. Only the identity 
(^|) is common to all models; the remaining bootstrap equations characterize the dilute 
Al model for different L = 3, 4, 6, 00. They have a purely geometrical origin in the 
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structure of the underlying algebra's root system and we conclude that just as the 
bound state structure of aifine Toda field theory can be encoded in Coxeter geometry, 
so can the string structure of the dilute A 3i 4 j6jQO models. The underlying Lie algebra 



also becomes manifest in another identity [|42| , f43| which can be built up using the more 
"fundamental" bootstrap identities ([761): 



rankg 



/•;('/■ ' -rH ) r ; ( " - — ) = f| n (u) (2 " A) >' • (77) 



1=1 



Here A is the Cartan matrix of the particular Lie algebra in question. This identity 
directly reflects the Dynkin diagram structure. Like the more general bootstrap identi- 
ties, ([FJ]) it points towards a relation between strings belonging to different excitations 
which we will investigate in a forthcoming publication |4T| . 



5.1.2 Spectral gaps and affine Toda mass spectra 

The excitation spectrum (|65| ) explicitly states that in regime 2 there are at most 8,7,6 or 
4 continuous bands of eigenvalues of the transfer matrix for L = 3, 4, 6, 00, respectively. 
Each band belonging to one of the leading excitations is separated from the ground- 
state by a gap. We now calculate the associated critical amplitudes and corrections 
to the scaling behavior (TO for all conjectured fundamental excitations flBSp. It will be 
worthwhile to do this in the framework of Coxeter geometry in order to reveal the Lie 
algebraic structure in the scaling corrections and relate them to the affine Toda mass 
spectra. 

As before we evaluate the excitations at the isotropic point u = 3 A/2 where the 
eigenvalues are real and set the parameter a in fl65|) to the value «o = (h + 2)t/tt giving 
the lowest eigenvalue in each band. We then find 

h # 4 (tl + ™ ; Id^r) 

Since we are interested in the behaviour of the correlation length in the critical limit, 
it is favorable to rewrite the above expression in terms of a power series in the elliptic 
nome using a standard identity for the logarithm of Jacobi's theta functions, 



= 4 ± ^-^^ t W«) - . (79) 



n=l 



From the coefficients in this expansion, which contain the Coxeter integers, we need 



now to extract the affine Toda masses. This has been done in the literature for 
L = 3,4, 6 case-by-case using trigonometric identities. It is worthwhile to review these 
calculations using Coxeter geometry to reveal the underlying Lie algebraic structure. To 



23 



this end we note the following formula for the residue of the inverse g-deformed Cartan 
matrix: 

(— l) s i v-^ / \ 7ra iP s , s 

= Sr~ E 2 ^ a sm IT 3 = ~T~^ ■ 80 

a=l « 

Here l<s</i — lisan exponent of the Lie algebra in question and we have used the 
identity fl67|). The second equality in (|80| ) involves the matrix 

P s = ^~ j d(((-2 + A)- 1 (81) 

|C-2cos^|= £ 

which is the orthogonal projector onto the eigenspaces of the (non-deformed) Cartan 
matrix, 

•7T C 

(P s ) 2 = p s and AP" = 4 sin 2 — P s . (82) 

The matrix elements of the orthogonal projector consist of the components of the nor- 
malized eigenvectors. Setting s = 1 we obtain the components of the Perron- Frobenius 
eigenvector which yields the (classical) affine Toda masses [[Tl[] 

P s k f l = m kmi . (83) 

Since all eigenvalues of the Cartan matrix are given in terms of the exponents, we 
conclude that the residue ( |80| ) vanishes if s is an integer which does not belong to the 
set of exponents. The final result for the correlation length and the leading term in the 
critical limit p — > therefore reads 

j. / \ i v^2/isin^r ™ \p\^ B m+rrik . . . 

Up)' 1 = E y— ti ' ,am § ~ 2/ > — ^ + • (84) 

- s sm - 1 - \ p \ h sin ^ 

Here s runs over the affine exponents, i.e. the finite exponents s modulo multiples of 
the Coxeter number, s = smodh. They are given in Table 5. 



Algebra 


exponents s 


E$ 


1,7,11,13,17,19,23,29 


E 7 


1,5,7,9,11,13,17 


E 6 


1,4,5,7,8,11 


D 4 


1,5,7,11 



Table 5. The exponents of particular finite Lie algebras, which determine the eigenvalues of 
the corresponding Cartan matrix. 



Res [A]~l ex 
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Note that only the odd exponents contribute to the sum in (I53J) and that we have chosen 
a normalization for the lightest mass different from mi = 1. Together with the result 
for the free energy density (|2TD the universal quantity (|3]) is expressed as 



Q ± (0)S ± (0y 



cos • 



2hm±mi 



15) 



The higher powers of the elliptic nome in the scaling functions (|84]) originate from ir- 
relevant operators of the minimal CFT Ml^l+i and the specific choice of the scaling 
variables. For the Ising model in a magnetic field, both contributions to the renormal- 



ization behavior have been analyzed in |0J . We therefore briefly comment on the scaling 
corrections in this case. 



The scaling functions for L = 3. Setting L = 3 in ( jS4]) we explicitly write out the 
first scaling corrections to the spectral gap of the dilute A 3 model 



= s^o)- 1 (i + \ P fl + \pfl + sM n + 



16) 



In comparison with |44| we note the absence of the analytic terms. This absence might 
be explained by a different choice of scaling variables, i.e. the identification of the elliptic 
nome with the magnetic field in the Ising model differs from |fH |. One has to keep in 
mind that the dilute Al model is not necessarily equivalent to the Ising model at the 
critical point but only resides in the same universality class. In contrast, the terms 



depending on the irrelevant operators of the identity family in M 3: 4 are present [44 |. 
The latter are related to the components TT,T 2 ,T 2 of the energy momentum tensor. 

From the expansion fl2~0|) we may also compare the scaling corrections for the singular 
part of the free energy density with those of [fHfl, 



\p\~™fs = Q ± (0) (l + Qi\p\% + Q 2 \pfl + Q 3 \pfl + -.) 

Once more we find all terms originating from the irrelevant operators belonging to the 
class of the identity field reported in jP[ . 



6 Conclusions 

In this article we have put forward the leading excitations of the dilute Al model in 
regime 2 as consisting of a hole and a two-string. Our proposal is for special values 
of L supported by numerical computations and further support comes from revealing 
the expected Lie algebraic structures £^,7,6 an d -D4 for L = 3,4,6 and L = 00, i.e. the 
unrestricted model. A crucial role in exhibiting these structures is played by analyticity 
arguments and Coxeter geometry which have allowed us to relate the zeroes in the exci- 
tation spectrum induced by the Bethe roots directly to the underlying Lie algebra and to 
cast the excitation spectrum of all four cases into a single formula (]65|) . This procedure 
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is more direct than using the non-linear integral equations of the thermodynamic Bethe 
ansatz, where one relates the corresponding integral kernels to the g-deformed Cartan 



matrix in Fourier space as has been done for L = 3 [23]. We note that the excitation 
spectrum ( |65D for Eg in fact matches the thermodynamic Bethe ansatz outcome of |23| 
and the result of the exact perturbation theory approach Q but we leave a closer in- 
vestigation of the proposed string structure for the higher excitations to future work 
|40|] . We believe that the one-to-one correspondence between the functional equations 
(f76|) and the bootstrap identities of the respective scattering matrices will be of central 
significance in gaining a thorough understanding of the allowed strings and a possible 
classification of the excitations. 

Our result for the spectral gaps and the critical amplitudes connected with the 
correlation length can now be compared against the result (^) from the thermodynamic 
Bethe ansatz analysis of the relativistic scattering matrices. We briefly report the result 
for the cases L = 3, 4, 6 where the associated scattering matrices can be written in the 



following universal formula (see e.g. |42|, |43j for a derivation): 



S(P) = exp r — 4 cosh t ([A]-K) sinh . (87) 
Jot V q / ii in 

We note that the g-deformed Cartan matrix [A] q which we frequently used in our analysis 
of the transfer matrix excitation spectrum also appears here. A similar calculation to 
that we have performed in section 4 in the context of elliptic functions, yields upon 
invoking formula (|4]) 

Q ± (0)5 ± (0) 2 - tan ^ 



\hm{ 



which is found to coincide with our result (|85|). Here m\ denotes the smallest component 
of the Perron-Frobenius eigenvector of the Cartan matrix, which can be identified with 



the classical mass. For general height values L our result ( j60| ) confirms the result in [T7 
based on the thermodynamic Bethe ansatz analysis of Smirnov's scattering matrices 
| III when the hole excitation is identified with the presence of a kink state. As the 



universal amplitude enters directly into the bulk vacuum expectation value of the energy- 
momentum tensor needed in the form factor approach to construct correlation functions, 
this underlines the intimate relationship between integrable lattice models and field 
theories. 

Another field theoretic aspect which needs to be mentioned is the interpretation of 
the particle spectra in light of the conformal structure at the critical point. For the 
special cases L = 3,4,6 where an underlying exceptional Lie algebra is present, the 
particle spectra manifest themselves in Roger s-Ramanuj an identities for the conformal 
characters ^7|, |5B |. Whether similar identities can be found for general L and 



further information on the conformal field theory can be extracted from the excitation 
spectrum are open challenges. 

We close by pointing out that further universal amplitude ratios can be derived from 
the critical amplitudes we have computed in this article. For example, provided that L 
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is even the elliptic nome can be identified with the reduced temperature up to a metric 
factor, p = — t = 1 — T/T c . This allows one to compute the following universal quantity 
involving the (singular) specific heat C s 0, 

R f = A^s ± (0), c s = -^~-r ft . (89) 

s ot z a 

Employing the concept of hyperscaling, the critical exponent a controlling the diver- 
gence of the specific heat near the critical point can be related to the critical exponent 
controlling the divergence of the correlation length: 

2z/ = 2-a = (l-A li2 )- 1 = i^^, Leven. (90) 



The coefficient A in fl89|) is easily obtained from the expansion of the singular free energy 
density fl2T|) to be 

A = 4:V3a(l-a)(2-a) ^~rrr ■ (91) 

sin f LT2 

The important observation in this context is that for L even the critical amplitudes S ± (0) 
of the correlation length above (p < 0) and below (p > 0) the critical temperature do 
not necessarily coincide. This is indeed the case for L — 4, i.e. the leading thermal 
perturbation of the tricritical Ising model, where the particle corresponding to the hole 
excitation is absent for p > and the heavier particle associated with the two-string 
gives the leading contribution |22[ . The different universal amplitudes above and below 
the critical temperature are related through the ratio fl62| ) which for L = 4 coincides 
with the previously reported result POL p2l [O] 



S + (0)/S-(0) = 1/2 cos § . (92) 
The absence of the hole excitation in depending on the sign of the elliptic nome should 



follow from symmetry arguments similar to those in the field theoretic context |49], 
|55| . To see this directly from the Bethe ansatz is an interesting problem for future 
investigations and requires a deeper analysis of the Bethe equations of the dilute Al 
model in the thermodynamic limit and of the string structure of the excited states. 
Finally, we expect that the analysis of the excitation spectrum via functional equations 
as outlined in this work can also be applied to the other regimes of the dilute A L model, 
though in a possibly modified form. The relevant universal amplitude ratios (Q) and 
(|39|) for regime 1, corresponding to the 2 i perturbation, have recently been reported 
in 1501. 
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A The eigenspectrum after the conjugate modulus 
transformation 



We briefly sketch how to prove the asymptotic behavior (|23"D in the ordered limit \p\ — > 1. 
Exploiting the conjugate modulus representation @ of Jacobi's theta functions the 
eigenvalue spectrum (|IDD can be rewritten as 



-e 2u ^wf) N A{w) =u 



E(x 4s /w)E(x 6s /w) 



-lN N 



E(x 4s )E(x 6 



ih 



-*s.E(x 2s w/wj) 
E(x 2s Wj/w) 



+ 



x 2s E(w)E(x 6s /w) 
wE(x 4s )E(x 6s ) 



N N 

n 



WjE(w/wj)E(x 6s Wj/w) 



LA E(x 2s Wj/w)E(x /is Wj/w) 



+ UJ- 



wE(w)E(x 2s /w) 
x 6s - 2r E(x is )E(x 6s ) 



iN N 



n 

i=i 



E(x 



8s-2r 



W : 



Wj/w) 



E(x 4s Wj/w) 



Here we have introduced the variables w = e ^~ u , s = tz§ and x = with r = 4t^| 
in regime 2 with p > 0. The no me q 2 = x 2r is suppressed in the notation. (Note that 
the parameters r and s and hence the variable x have been rescaled from those used 
by the authors of |57|] in their various papers.) Assuming that for the groundstate all 
Bethe roots lie on the unit circle we derive for < r < 1 (i < 1), 



e -r A (w) 



^o n ^ 

n 



n ■ x+y < M < A 

A < u < 2A 



„2sV 



U X— U W n r + (rfer II Wn + 00 Q ]J W r l 



-1 



2A < u < 3A 



and 



"AoM 



3jS * 

-W r 



N 



+ 



>x 2s \N 



, 0< M<Z ^. 



We therefore have the following asymptotic behaviour as A" — ► oo: 



2n 2 iV . , N 

e ^ A (w) 



u; n w « 



uj 1 n ^ 



< Reu < A 

A < Reu < 2A 
2A < Re u < 3A 
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We thus conclude that in the thermodynamic limit iV ^> 1, 



Ao(«) 



LOq 



i9i(2A~m)^i(3A-m) 
j?i(2A)#i(3A) 

i?i(m)i?i(3A-m) 
i?i(2A)i?i(3A) 



iV 



} G (u)[l + o(e- N )] 
[1 + °(e-")] , 



< u < A 
A < it < 2A 

^o 1 !^®^} G ( M -3A)- 1 [l + o(e^)], 2A < u < 3A 



Go(m-A) 
G (n-2A) 



Since the excited states A(w) differ only by a finite number of Bethe roots with non- 
vanishing real part we obtain analogous asymptotic behaviour in the thermodynamic 
limit for them also. This proves (|2"B|) and therefore the functional equation 



B Coxeter geometry 

Let g denote in the following a simple Lie algebra and denote by {«i, ...,a n } a set of 
simple roots which span its root system, i.e. the eigenvalues of the Cartan subalgebra 
in the adjoint representation. Each simple root on can be interpreted as element in 
an Euclidean vector space R n of dimension n =rankg. It naturally defines a reflection 
R n — > M n at its associated hyperplane by setting 

v^Knv:=v- 2 ^' "V on , veR n . (93) 

The transformations <7j are called simple Weyl reflections and generate the Weyl group 
W. Their action on simple roots is described in terms of the Cartan matrix associated 
with q, 

2(cti } aj} 

y Gi j - Gt j J 

There exists a longest element in W in the sense that it is built up from a maximal 
number of simple Weyl reflections. It is called the Coxeter element or transformation 
and is defined by the product over all simple Weyl reflections, o = a±a 2 ■ ■ ■ cr n . Clearly 
this definition depends on the ordering of the reflections, which is a matter of choice. 
The Coxeter element is therefore only defined up to conjugacy within the Weyl group. 
Following [[31] a unique Coxeter element can be determined by introducing the concept 
of bicolouration for Dynkin diagrams. To every vertex in the Dynkin diagram we assign 
a colour Cj = ±1 such that two vertices linked to each other are differently coloured. 
This bicolouration polarizes the index set A = {1, n} into two subsets A± and allows 
unambiguous specification of the following Coxeter element: 

a := <7-<J + , a± := ]~J t7» . (94) 



Any Coxeter element shares the following properties 51 : 
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(CI) The Coxeter element fixes no non-zero vector. 



(C2) It is of finite order, a h = 1, where h is the Coxeter number defined as the sum 
over the Kac or Dynkin labels, 

rank g 

h := 1 + n i ■ ( 95 ) 

i=l 

Thus, the Coxeter element permutes the roots in orbits of length h. 
(C3) The eigenvalues of a are of the form 

exp-7-, j = 1, rankfl 

where the characteristic set of integers 1 = sx < s 2 < ... < s n = h — 1 are called 
the exponents of the Lie algebra g satisfying the relation s n+ i_j = h — Sj. 

The action of the Coxeter element on the root system splits the latter into orbits. 
A convenient choice of representatives of these orbits are the ,, coloured" simple roots, 
7j = CjO!j. The Coxeter orbits fli defined as 

■= { a x li :l<x<h} (96) 

satisfy the crucial property that they do not intersect, i.e. f^nfij = 0, and are exhaustive 
on the set of roots. Moreover, all 7/s lie in different orbits and all elements in one orbit 
are linearly independent f3]| . Thus, the coloured simple roots constitute a complete set 
of representatives for the Coxeter orbits fli. 
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